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Large-scale magnetic fields in the inflationary universe 
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The generation of large-scale magnetic fields is studied in inflationary cosmology. We consider 
the violation of the conformal invariance of the Maxwell field by dilatonic as well as non-minimal 
gravitational couplings. We derive a general formula for the spectrum of large-scale magnetic fields 
for a general form of the coupling term and the formula for the spectral index. The result tells 
us clearly the (necessary) condition for the generation of magnetic fields with sufficiently large 
amplitude. 



I. INTRODUCTION 
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' It is well known that there exist magnetic fields with the field strength ~ lO^^G on 1 — lOkpc scale in galaxies of all 
^ ! types (for detailed reviews of cosmic magnetic fields, see [U, H, H, 0, S H, 0| ) and in galaxies at cosmological distances 
p^J ' H- Furthermore, magnetic fields in clusters of galaxies with the field strength 10~^ — 10~^G on lOkpc— IMpc scale 
CnI ] have been observed Q . It is very interesting that magnetic fields in clusters of galaxies are as strong as galactic ones 
and that the coherence scale may be as large as ~Mpc. The origin of these magnetic fields is not well understood yet. 
Galactic dynamo mechanisms [lO[ have been proposed to amplify very weak seed magnetic fields up to ~ 10~^G. In 
fact, however, galactic dynamo mechanisms require initial seed magnetic fields to feed on, and the effectiveness of the 
' dynamo amplification mechanism in galaxies at high redshifts or clusters of galaxies is not well established. 
, Proposed scenarios for the origin of cosmic magnetic fields fall into two broad categories. One is astrophysical 
^— I ■ processes, e.g., the Biermann battery mechanism [ll| and the Weibel instability which is a kind of plasma 
T— I ' instabilities, and the other is cosmological processes in the early Universe, e.g., the first-order cosmological electroweak 
phase transition (EWPT) [l^, quark-hadron phase transition (QCDPT) (see also and the generation of 

the magnetic fields from primordial density perturbations before the epoch of recombination piil ITtI . fisl . [l9| . It is 
(~| : difhcult, however, for these processes to generate the magnetic fields on megaparsec scales with the sufficient field 
i' strength to account for the observed magnetic fields in galaxies and clusters of galaxies without requiring any dynamo 
Q \ amplification. 

^ ' The problem is essentially the same as the problem of the large-scale structure formation for which one needed the 
' adiabatic curvature perturbation of non-negligible amplitude on superhorizon scales, the so-called horizon problem. 
d • Then, just as this problem is solved by the inflationary cosmology, it is natural to look for the possibility of generating 
'>-^ ' such a large-scale magnetic field during inflation [20]. However, the Friedmann-Robertson- Walker (FRW) universe is 
conformal fiat and the Maxwell theory is conformal invariant. Therefore, the conformal invariance must be broken 
to generate non-trivial magnetic fields.^ Various conformal symmetry breaking mechanisms have been proposed so 
^ ; far [20, 22, 23, 24, 25, 26, 27, 28, 29, HO, iH, [32, 33, 34]. 

Bamba and Yokoyama studied the case of the dilaton electromagnetism in Ref. [35| . where the Lagrangian is in 
the form /($)F^. The existence of the dilaton is motivated by string theory, and such a coupling is reasonable in 
the light of higher-dimensional theories. The case of the coupling function /($) = e"^''*, where k = VSttG, was 
first analyzed by Ratra ^2^, in which the inflaton and the dilaton are identified. Bamba and Yokoyama then showed 
that magnetic fields with the current strength as large as 10~^''G on 1 Mpc scale could be generated [1^. For this 
to be the case, however, they had to introduce a huge hierarchy between the coupling constant of the dilaton to the 
electromagnetic field A and the coupling constant A in the dilaton potential, A/A 3> 1, in order for the spectrum of 
the generated magnetic fields to be sufficiently red so that the amplitude on large scales can be sufficiently large. 

As a possible solution to the above huge hierarchy between A and A, Bamba and Yokoyama proposed a new 
scenario [s^ by taking account of the effects of the stringy spacetime uncertainty relation (SSUR) [s^l- They found 
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^ In Ref. [2ll . the effect of the breaking of the conformal flatness due to scalar metric perturbations at the end of inflation is discussed. 
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that the SSUR of the metric perturbations could lead to magnetic fields with a nearly scale-invariant spectrum even 
if A and A are of the same order of magnitude, and the amplitude of the generated magnetic fields could be as large 
as 10^^"G on IMpc scale at the present time. 

As another mechanism to break the conformal invariance, non-minimal gravitational couplings were proposed by 
Turner and Widrow [20] . Many of the couplings they introduced break not only the conformal invariance but also the 
gauge invariance. The one that looks theoretically more appealing is of the form RF^ , where R is the Ricci scalar. 
Such a term is known to arise in curved spacetime due to one-loop vacuum-polarization effects ^3§\ . 

In the present paper, we consider the case when there are not only dilatonic and curvature couplings but also any 
possible coupling of the form /(<&, R, -- where the dots (• • • ) can be anything. This paper is organized as follows. 
In Sec. II, we describe our model and derive the equations of motion. In Sec. Ill, we consider the evolution of the U{1) 
gauge field and that of the magnetic fields in a generic slow-roll inflation model. We derive formulas for the spectrum 
of the generated magnetic field and the specral index expressed very generally in terms of the coupling function 
/($, i?, • • • ) without specifying its explict form. These general formulas tell us clearly the necessary condition for the 
generation of a sufficiently large amplitude magnetic field on large scales. Finally, Sec. IV is devoted to conclusion. 



II. MODEL 



We consider a conformal violationg Maxwell theory with the action, 

d^xy^ i?, • • • ) F.^Ff^--^ , (2.1) 

where / is an arbitrary function of non-trivial background fields at the moment. Naturally the dots (• • ■ ) in the 
argument of / will include the inflaton if the dilaton $ does not play the role of the inflaton. 
From the action (|2.ip . the electromagnetic field equation is given by 

^ -.d^ [V^IF^^''] - 0. (2.2) 



V-9 

We assume a spatially flat universe with the metric 

ds^ = -df + a^{t)dx^ = a'^{'q){~drf + dx^), (2.3) 

where a is the scale factor and 77 is the conformal time. We consider the evolution of the U{1) gauge field in this 
background. Its equation of motion in the Coulomb gauge, Ao(t, a;) = and d^Aj{t, x) = 0, becomes 

f i\ ■ 1 (3) 

A,{t,x)+ \H+j\A,{t,x)^ —AA,{t,x) = 0, (2.4) 
where H — d/a, and a dot denotes a time derivative, ' ~ d/dt. 



III. GENERIC SLOW-ROLL INFLATION 



In this section, we consider the evolution of the U{1) gauge field and that of the magnetic fields in generic slow-roll 
inflation. 



A. Spectrum from quantum fluctuations 

To begin with, we shall quantize the U{1) gauge field A^(i, cc). It follows from the electromagnetic part of our 
model Lagrangian in Eq. p.ip that the canonical momenta conjugate to A^{t,x) are given by 

TTo = 0, TT, = Ia{t)A,{t,x). (3.1) 
We impose the canonical commutation relation between Ai{t, x) and TTj{t, x), 

[A{t,xln,it,y)]=^ J -g^e^f— > - ^) , (3.2) 
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where k is comoving wave number and k = \k\. From this relation, we obtain the expression for Ai(t,x) as 

— y \b{k,a)e,ik,a)Ait,k)e"'- +b\k,a)eUk,a)A*{t,k)e-'''- 



(3.3) 



where ti(k^a) (a — 1,2) are the two orthonormal transverse polarization vector, and b{k,a) and 6^(A;,(t) are the 
annihilation and creation operators which satisfy 



b{k,a),P{k',a')] ^ 5a,a'6^{k ~ k'), [^(A;, ct), 6(A;', ct') 1 ^ \b\k,a),b\k\a') 



= 0. 



It follows from Eq. 



that the mode function A{k, t) satisfies the equation 
A{k, t)+\H+- \ A{k, t) + ^ A(fc, t) = 0, 



(3.4) 



(3.5) 



with the normalization condition. 



A{k,t)A*{k,t) - A(k,t)A*(k,t) = — , 

la 



Replacing the independent variable t to 77, we find that Eq. 



becomes 



A 



r 



"{k,v) + jA'{k,7^) + k^A{k,v)^0, 



(3.6) 



(3.7) 



where the prime denotes differentiation with respect to the conformal time ij. 

Although it is impossible to obtain the exact solution of Eq. (|3.7p for the case when / is given by a general 
function of 77, we can obtain an approximate solution with sufficient accuracy by using the WKB approximation on 
subhorizon scales and the long wavelength approximation on supcrhorizon scales, and match these solutions at the 
horizon crossing. 

In the exact de Sitter background, we have a = \/{—Hrj) where H is the de Sitter Hubble parameter. Hence the 
horizon-crossing, which is defined by _ff = fc/a, is given by ^krj = 1. The subhorizon (supcrhorizon) scale corresponds 
to the region k\ri\ ^ 1 (fc|?7| <C 1). We expect this to be also a sufficiently good definition for the horizon crossing for 
general slow-roll, almost exponential inflation. 

Then the WKB subhorizon solution is given by 



Ain{k,ri) 



1 



where we have assumed that the vacuum in the short-wavelength limit is the standard Minkowski vacuum. 
The solution on supcrhorizon scales can be obtained by using the longwavelength expansion, 

^out =Ao(?7)-ffc2Ai(r,)+0(fc4). 
Substituting Eq. ((3J| into Eq. ((37)) . we obtain 



(3.8) 



(3.9) 



K + jK = 0, 



A'l + y a; + Ao 



0. 



(3.10) 
(3.11) 



Let the two independent solutions for Aout be u and v. For definiteness, we set the boundary conditions to be m 1 
and z; ^ as 77 ^ T^R, where rj^ is the conformal time at the time of reheating after infiation. From Eq. p.lOp . the 
lowest order solutions are given by 



Mo 



1, 



■zzdfj . 



(3.12) 
(3.13) 
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The solutions accurate to 0{k'^) can be obtained by the iteration using Eq. (|3.1ip . We find 



uq + k-UQ 
vo + 



" drj" 



With these solutions, the general superhorizon solution is given by 

^out = Cu + Dv, 



lin") ' 

" drj" 



(3.14) 



(3.15) 



where C and D are constant. 

We match the subhorizon solution (|3.8p with the superhorizon solution (I3.15P at the horizon crossing t] ^ rjk ~ —/k 
to determine the constants C and D. The junction conditions are 



Ain iVk) = Aout iVk) 



(3.16) 
(3.17) 

The determination of C and D accurate to 0{k^) are given in Appendix. It turns out that O^k"^) corrections are small 
unless the function I is non-monotonic and sharply decreases on superhorizon scales. Hence we neglect the 0{k^) 
corrections here and approximate Aout by the lowest order solution, 



Aout = C{k) + D{k) 



nn 



1 



Substituting Eqs. ([3^ and ([STS]) into Eqs. ([3T6l) and (|3T7| . we obtain 



C{k) = 



-1/2 



2k 



iki 



r 1 ^ ' 



-ikrj 



(3.18) 



(3.19) 



D{k) 



(3.20) 



V=Vk 



Neglecting the decaying mode solution, it then follows from Eqs. (|3.18|) that |^(fc,77)p at late times is given by 



\A{k,^)\' = |C(fc)p = 



1 



2kl{rjk 

The proper magnetic field is given by 



(3.21) 



(3.22) 



where Bi{t,x) is the comoving magnetic field, and e.yfc is the totally antisymmetric tensor (ei23 = 1)- Thus the 
spectrum of the magnetic field is given by 



|i3p™p-(fc,r;)p = 2^|A(fc,r;)p = 2^|C(fc)p 



(3.23) 



where the factor 2 comes from the two polarization degrees of freedom. Thus the energy density of the generated 
magnetic field per unit logarithmic interval of k is given by 

1 47rfc3 



PBik,ri) 



2 (27r): 



|S™(fc,ry)p/(ry) = ffl^^/(r;), 



and the density parameter per unit logarithmic interval of k after reheating becomes 

„ . PBjk.m) m l5k\Cik)l ( , M 

ilB(k,ri) = — — — — r = ^4 4 7—I[V), PimJ^N^si^ 



30 



R > 



(3.24) 



(3.25) 



where Tr is the reheating temperature, or is the scale factor at reheating, and Ncs is the effective massless degrees 
of freedom (2 for photons) which are thermalized at reheating. The spectral index of Qsik) is then given by 



_ dlnneik) _ ^ 
dink 



dlnfc|C(fc)p 
din fc 



(3.26) 
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B. Semi-quantitative estimate 

Using Eq. (|3.21[) for |C(fc)p, one may derive a general expression for the density parameter and the spectral index of 
the magnetic field in terms of /. However, the result will be very complicated and hence will not be very illuminating. 
To gain more insight into semi-quantitative nature of the generated magnetic field, let us therefore a specific form for 
the function /. We set 

(3.27) 

where rj^ is some fiducial time during inflation, and a is a constant. As we shall see shortly, in order to obtain a 
cosmologically interesting result, we need a monotonically increasing / during the stage of inflation when the relevant 
comoving scales leave the horizon. So, we assume a > in the following. Then |Cp in Eq. p.2ip can be explicitly 
evaluated to be 



k\C\'= ^ 



21{vk) 



_ a + 2i 
2{a + l) 



2/(?7fe) ' 



(3.28) 



where kr]k = 0(1) and C is a constant of order unity. 

Then, inserting the above equation to Eqs. (|3.25p and p.26p . we obtain the explicit expressions for the density 
parameter and the spectral index at present rjo'- 

f^s(fc,r;o)--4-^— -— ?^ , nB = 4-«, (3.29) 



where 1(7]^) cx k" and /(?/o) = 1 have been used. We see that for a spectrum with a sufficiently large amplitude, 
must be sufficiently small, and for a spectral index close to a scale-invariant one, we need a ~ 4. 

To proceed further to make a more quantitative estimate, let us assume that the reheating takes place almost 
instantaneously so that rjn can be identified with the conformal time at the end of inflation, and /(ry) is given by 
Eq. ((3?27)) until the end of inflation. Noting that 0^77^ « H^'^ and 3ff| = Pr/M^j, where Mpi = l/VSnG, from 
Eq. ((3?25ll we obtain 

We see that higher the reheating temperature is, larger the density parameter of the magnetic field becomes. But given 
the fact that Tr cannot be made arbitrarily large, there is a limit to this mechanism. Instead the most important 
point of the above formula is the presence of the factor 1//(77r). It means one can make fls arbitrarily large by 
choosing /(77R,) arbitrarily small. An explicit example of this case was discussed in [35| . where / was a function of a 
dilaton which is assumed to be rolling even after the end of inflation. It was found that the present magnetic field 
with rms amplitude as large as 10~^ Gauss can be obtained. 

Of course, one should be aware of the fact that it is extremely difficult to construct a realistic model that gives 
both a rapidly increasing coupling function / and a very small value of / at reheating at the same time. Even just 
one of these two conditions is very hard to realize. 

Before closing this section, we mention one important fact. From the above result, one might have an impression 
that there is an extra enhancement of the magnetic field from the post-inflationary stage. If this were true, it would 
mean the existence of a new enhancement mechanism. But it is not so. The reality is that the whole enhancement 
factor is entirely determined by the amplitude of / at horizon-crossing during inflation. This can be seen from the 
form of fis given in Eqs. (|3.29p in its original form, 

Tga^ 2A'eff7r4 I(rik) 

The above expression clearly shows that it is indeed the value of / at horizon-crossing that determines the final 
amplitude as well as the spectral index. 



IV. CONCLUSION 



In the present paper we have studied the generation of large-scale magnetic fields in infiationary cosmology, breaking 
the conformal invariance of the electromagnetic fields by introducing a coupling IF^ with / depending on non-trivial 
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background fields that vary in time. We have derived a general formula for the spectrum of large-scale magnetic fields 
for a general form of the coupling term and the formula for the spectral index. The result shows the necessary condition 
for the generation of magnetic fields with large amplitude. Namely, the coupling function / must be extremely small 
and rapidly increasing in time during inflation. Such a coupling seems very unnatural at first sight. Nevertheless, 
considering the recently much discussed issue of the string landscape, which tells us that infinitely many different low 
energy theories are possible |3i, the naturalness may not be crudal, and it may worth searchfng for a model that 
realizes a successful inflationary universe and satisfies the necessary condition for the generation of a large magnetic 
field. 
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APPENDIX 



In this appendix, we determine the constants C and D in Eq. (j3.15[) by taking into account the terms up to order 
Substituting Eq. ([STTS]) into Eqs. ([3?T6l) and (pTZ]) . wc find 



Am 

AL 



M 



n=Vk 



where 



M 



h (1 - fc^/i 



Ph/I [-l + k\h+hh)\/I 
Here, /i(?7), l2{'>]), and Isi^r]) are defined by 



hiri) ^ 



I in) 



dfj, 



hirj) EE / Iif^)dfj. 



From Eqs. (|XT|) and (|X2|) . we obtain 



= M~ 



1 ( Ain 

AL 



where 



1 f l-e{h+ hh) Ih (1 - fc'/i) 



(A.l) 
(A.2) 

(A.3) 
(A.4) 
(A.5) 

(A.6) 

(A.7) 
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